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ABSTRACT 

Observations by the Large Area Telescope (LAT) on the Fermi mission of diffuse 
7-rays in a mid-latitude region in the third quadrant (Galaetic longitude / from 200° 
to 260° and latitude |6| from 22° to 60°) are reported. The region contains no known 
large molecular cloud and most of the atomic hydrogen is within 1 kpc of the solar 
system. The contributions of 7-ray point sources and inverse Compton scattering are 
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estimated and subtracted. The residual 7-ray intensity exhibits a hnear correlation 
with the atomic gas column density in energy from 100 MeV to 10 GeV. The measured 
integrated 7-ray emissivity is (1.63 it 0.05) x 10~'^^ photons sr~^ H-atom~^ and 
(0.66 ±0.02) X 10-26 photons s'^ sv~^ H-atom^i above 100 MeV and above 300 MeV, 
respectively, with additional systematic error of ~ 10 %. The differential emissivity from 
100 MeV to 10 GeV agrees with calculations based on cosmic ray spectra consistent with 
those directly measured, at the 10 % level. The results obtained indicate that cosmic 
ray nuclei spectra within 1 kpc from the solar system in regions studied are close to the 
local interstellar spectra inferred from direct measurements at the Earth within ~ 10 %. 

34 Subject headings: cosmic rays - diffuse radiation - gamma rays: observations 



1. Introduction 



36 The diffuse high energy 7-ray emission (E > 30 MeV) has been interpreted to be a superposi- 

37 tion of 7-rays produced via interactions between cosmic rays (CRs) and interstellar matter, inverse 
3B Compton (IC) scattering of interstellar soft photons off CR electrons, and the extragalactic diffuse 

39 7-ray emission. The first component, if distinguished from the others, will enable using high-energy 

40 7-ray observations for the study of the distribution of CRs and the interstellar medium. The distri- 

41 bution of neutral atomic hydrogen (H I) is traced by 21 cm line surveys and the molecular hydrogen 

42 distribution is derived indirectly using 2.6 mm line observations of carbon monoxide (CO). The total 

43 gas column density can also be traced indirectly from extinction and reddening by dust. Thus the 

44 spectrum and the flux of CRs can be obtained from sufficiently sensitive observations of high energy 

45 7-rays. The observation of diffuse 7-rays away from the Galactic plane (Galactic latitude |6| > 10°) 

46 is suitable for studying local CRs, since diffuse 7-rays in such regions are less affected by contam- 

47 inati on from strong point sources, arid mo st of the gas along the line of si ght is local. The SA S-2 



(e.g.. iFichtel. Simpson. &: Thompson! 1 19781 ) and COS-B observations (e.g.. iLebrun et al.lll982l ) in- 



dicated a correlatio n between the 7-ray int ensiti es and the total gas column densities at medium 
Galactic latitudes. ISreekumar et al.l (jl998l ) and IStrong. Moskalenko. Sz Reimerj (j2004l ) showed a 
good correlation between the 7-ray intensities and model calculations in their analyses of the ex- 
tragalactic diffuse emission observed by EGRET onboard the Compton Gamma-Ray Observatory. 
Despite these early studies, the flux and spectrum of local CRs deduced from 7-ray intensity remain 
uncertain, due to the possible contamination from unresolved point sources and uncertainties in 
modeling the IC contribution because of the large scale height of CR electrons and the reprocessing 
of the interstellar radiation by dust. Although CR nuclei in the vicinity of the solar system are 
thought to have spectral distributions and inter isities similar to those measured at the Earth as 
reported by a number of 7-ray observations (e.g. iHunter et al.lll997l : iDigel et al.ll200ll ). data above 
1 GeV, which are crucial to distinguish CR nuclei spectra from that of CR electrons, have not been 
good enough due to the limited photon statistics and relatively limited energy coverage of these 
early missions. 
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The situation has improved significantly with the recent launch of the Fermi Gamma-ray Space 
Telescope on 2008 June 11. The Fermi LAT (Large Area Telescope) has a sensitivity that is more 
than an order of magnitude better than that of EGRET and enables resolving point sources and 
studying the diffuse 7 -rays with unprecedented se nsitivity. Recent advances of a CR propagation 



code GALPROP (e.g.. lStrong &: Moskalenkdll998l ). which had been developed through comparisons 



with the EGRET data, allow us to predict and subtract IC emission and correlate 7-ray emission 
with interstellar matter more accurately. 

In this paper, we present Fermi LAT observations of diffuse 7-rays in a mid-Galactic latitude 
region in the third quadrant (Galactic longitude / from 200° to 260° and |6| from 22° to 60°). As 
discussed in the following sections, most of the gas along the line of sight is local, nearby on the 
scale of the Milky Way. The contribution from IC emission is only about 10 % of the total diffuse 
emission and the LAT has already resolved five times as many 7-ray point sources as previous 
missions in this region. These facts enable us to evaluate the local CR flux and the spectrum with 
small systematic uncertainty. 



Observation and Data Reduction 



The LAT is the main instrument of the Fermi Gamma-ray Space Telescope. It consists of 
4x4 modules (towers) built with tungsten foils and silicon microstrip detectors to measure the 
arrival directions of incoming 7-rays, and a hodoscopic cesium iodide calorimeter to determine 
the photon energies. They are surrounded by 89 segmented plastic scintillators serving as an 
anticoincidence detector to reject charged particle events. Details of the LAT instrument and 
pre-launch expectations of the performance can be found in lAtwood et al.l (120091). T he excellent 



sensitivity of the LAT is exemplified by initial publications such as 



Abdoetal 



torn 



Routine science operations with the LAT began on 2008 August 4. We have accumulated events 
from 2008 August 4 to 2009 January 31 to study diffuse 7-rays. During this time interval the LAT 
was operated in sky survey mode nearly all of the time; in this observing mode the LAT scans 
the sky, obtaining complete sky coverage every two orbits and relatively uniform exposures over 
time. We used the standard LAT analysis software, ScienceTools version v9rll, and applied 
the following event selection criteria: (1) events have the highest pr obability of being photons, i.e., 
they are categorized as so-called diffuse class (lAtwood et al.ll2009l ). (2) the reconstructed zenith 
angles of the arrival direction of photons are selected to be less than 105°, in order to exclude 
periods where the Earth enters the LAT field of view, and (3) the center of the LAT field of view 
is within 39° from the zenith in order not to include the data taken in the pointed observation 
mode, because it has increased contamination from Earth albedo 7-rays. We also eliminated the 
period of time during which the LAT detected two bright GeV-emitting GRBs, i.e., GRB080916C 



available from the Fermi Science Support Center ( |http : //f enni . gsf c . nasa . gov/ ssc/ [ ) 
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(|Abdo et al.ll2009al ) and GRB081024B (|Omodei et al-lbooa ). We then generated count maps (using 
gtbin in ScienceTools) and exposure maps (using gtexpcube) in 13 logarithmically-sliced energy 
bins from 100 MeV to 9.05 GeV. A post-launch response function P6_V3_DIFFUSE, which was 
developed to account for the 7-ray detection inefficiencies that are correlated with trigger rate, 
was used in exposure calculations. These count and exposure maps we re prepared in Ca. rtesian 
coordinates in 0.5° x 0.5° binning, and then transformed into HEALPix[§ ( Gorski et al. 2005 ) equal 
area sky maps of order=7. They are used below to correlate the 7-ray intensities with the column 
densities of atomic gas along the line of sight; the 7-ray intensity is calculated as the ratio of the 
counts and the exposures for each energy bin. 



3. Data Analysis 

3.1. Subtraction of Inverse Compton and Point Sources 

To distinguish 7-rays produced in the interstellar medium from others, we referred to the GAL- 
PROP prediction of IC emission and an LAT source list for 6 month data|f|. This list w as produced 



109 using a similar procedure used to obtain the LAT Bright Source List described by lAbdo et al 



no (|2009bl ). It covers the same period of time as that of our data set and contains 740 point sources 

111 with significance more than 5 a. We adopted positions and spectral parameters from this list (single 

112 power-law model in 100 MeV-100 GeV) t o estimate and subtract the photons from point sources to 



113 diffuse 7-ray emission. GALPROP (e.g., IStrong &: Moskalenkdll998l ) is a set of programs to solve 



114 the CR transport equation within our Galaxy and predict the 7-ray emission produced via inter- 
ns actions of CRs with interstellar matter (nucleon-nucleon interaction and electron bremsstrahlung) 
116 and soft photons (IC scattering) . IC emission is calculated from the distribution of (propagated ) 



117 electrons and the interstellar radiation fields developed by Moskalenko. Porter. &: Strong (|2006l ). 
iiB Here we adopted the IC model map with version 54_5g;Xvarh7S B, w hich was used in a nother 



Fermi LAT paper to study the diffuse 7-ray emission in 10° < |6| < 20° (|Abdo et al.1 l2009dl ^ . The 



120 CR electron spectrum is adjusted to agree with the directly-measured pre-Fermi spectrum in this 

121 GALPROP model. In order to minimize the uncertainty of the contribution from IC emission on 

122 the diffuse 7-ray spectrum, we selected sky regions away from the Galactic center; the lower CR 

123 electron fluxes and interstellar radiation field will result in dimmer IC emission than that toward 

124 the Galactic center. We chose the third quadrant. Galactic longitude / from 200° to 260° and the 

125 Galactic latitude h from —60° to —22° and from 22° to 60°. The region is free of known large 

126 molecular clouds; Orion molecular clouds (Orion A and Orion B) and Monoceros molecular cloud 

127 complexes are located in the region I from 200° to 220° and b from —10° to —20°, and the Tau- 



^|http : //healpix . jpl . nasa . gov' 

^internally available to the LAT team 

4 



The GALPROP galdef ID of this version is available at the website jhttp: //galprop. Stanford. edu| 
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128 rus /Perseus molecular clouds are in I fro m 150° to 185° (e.g.. lDame. Hartmann. &: Thaddeusll200ll : 



Digel et al.lll999l : iDigel &: Grenieij 120011 ). Therefore the region described is suitable for correlating 



130 the 7-ray intensities with the local atomic-gas column densities. 

131 In Figure 1, we show 7-ray count maps above 100 MeV. There are 52 sources in the LAT 

132 6 month source l ist in our region of i nterest, more than five times as many sources in the third 

133 EGRET catalog (jHartman et al.lll999l ) in this region (nine sources). The diffuse 7-ray spectrum, 

134 after masking sources with circular regions of 1° radius, is shown in Figure 2. Atomic hydrogen 

135 column density maps of the same region (see Section 3.2 for details) are given in Figure 3. In 

136 Figure 2 and figures shown hereafter (Figures 4-6), the 7-ray intensities or CR fluxes multiplied by 

137 E'^ (where E is the center of each energy bin in logarithmic scale) are presented. Also presented in 
13B Figure 2 is the contribution from IC emission predicted by GALPROP, and the spillover from point 

139 sources outside the mask regions estimated (using gtmodel) by the spectral parameters given in 

140 the source list. Both the flux of estimated IC emission and the residual point source contribution are 

141 less than 15 % of the total diffuse emission above 100 MeV. We thus conclude that the uncertainty 

142 due to the IC and point source contributions is negligible after we subtract them from 7-ray data. 

143 Hereafter we analyze diffuse emission after masking point sources and subtracting IC emission and 

144 the residual contributions from point sources. 



3.2. Atomic Hydrogen Map 



146 Column densities A^(HI) of atomic hydrogen gas were calculated from existing radio surveys 

147 of the 21 cm line of HI. We used the Leiden/Argent ine/Bonn (LAB) Survey which merges the 



Leiden /Dwingelo o Survev (jHartmaiin &: Burton 



149 tronomia Survey (lArnal et al 



2OO1 



I 



Bajaja et al 



19971) with the Instituto Argentino de Radioas- 



2OO5I ) and covers the entire sky. We applied an 



150 optical depth correction under the assumption of a uniform spin t emperature of 125 K and the cos- 

151 mic microwave background intensity at 1420 MHz of 2.66 K (e.g.. iHunter et al.lll994l ). The derived 

152 H I column density ma ps of our region of interest are shown i n Figure 3. Although major CO sur- 

153 veys such as the one bv lDame. Hartmann. &: ThaddeusI (|200ll ) do not cover the region we analyzed, 

154 no large molecular cloud is known there and the molecular gas con tribution is expected to be small 

155 due to the moderately high Galactic latitude. See the discussion bv lDame. Hartmann. &: Thaddeus 

156 (|200ll ) for the completeness of their survey. Hereafter we assume that all the gas is in atomic form 



157 and traced by 21 cm radio surveys. Column densities of H I in our region range from 1 x 10 cm~ 

15B up to 18 X lO^'' cm~^ and the optical depth correction is rather small; the increase of the col- 

159 umn densities from those for the optically thin case (infinite spin temperature) is < 10 % in most 

160 directions. On the assumption of a Galactic rotation curve by IClemensI ()1985l ) for the case of 

161 Rq = 8.5 kpc and = 220 km s^^ (where Rq and are the Galactocentric radius and the orbital 

162 velocity of the local group of stars, respectively), we infer that, in almost every direction in our 

163 region, more than 80 % of the HI along the line of sight is within 1 kpc of the solar circle. Further - 

164 more, by referring to the vertical density distribution of HI given by iDickev &: LockmanI (|l990l ). 
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165 we can conclude that more than 85 % of atomic gas in the hne of sight is within 1 kpc of the solar 

166 system for |6| > 22°. 



3.3. Correlation of 7-ray Intensities and Gas Column Densities 



The LAT point-spread function (PSF) strongly depends on the photon energy (e.g. JAtwood et al 



2OO9I ) and the energy dependence of the angular size needs to be taken into account in data analy- 



sis. We convolved t he map of H I column densities obtained as described in Section 3.2 using the 



GaDGET package (jAckermann et al.ll2008l ) with the LAT PSF for each of our energy bins; in the 
convolution we used the all sky map to take account of the contribution from outside the region for 
the analysis. Since the typical angular size of the variation of column densities is a few degrees (see 
Figure 3), only maps for the lowest energy bands (less than a few hundred MeV) are noticeably 
smeared. 

The 7-ray intensities, after masking point sources with 1° circular regions and subtracting the 
IC emission and the residual point source contributions, are correlated with the H I column densities 
in each energy band. Both the 7-ray intensity map and the H I column density map were prepared 
in HEALPix equal area sky maps of order=7, whose pixel size in solid angle is 6.39 x 10~^ steradian 
and is close to that of 0.5 x 0.5 deg^. We found a linear relationship between A^(HI) and residual 
7-ray intensities for energies from 100 MeV to 10 GeV. Above 10 GeV the correlation is limited 
by photon statistics. Figure 4 shows the correlation between 7-ray intensities and the HI column 
densities for four representative energy bands. The linear correlation indicates that point source 
contributions are successfully subtracted and residual 7-rays mostly originate from interstellar 
atomic gas through interactions with CRs, plus isotropic diffuse component (extragalactic diffuse 
7-rays, the residual particle background, and a possible residual of IC emission). 

By fitting the correlation in each energy band with a linear function using a minimization, 
we obtained the intensity of the isotropic diffuse component and the emissivity of atomic gas as 
the offset and the slope, respectively, as summarized in Table 1. Making the mask region larger, to 
3° radius, gives consistent fit parameters within statistical errors, confirming that the contribution 
of point sources is well modeled and subtracted. The obtained isotropic diffuse comp onent (offset 



192 column in Table 1) agrees within 10-20 % with the "Isotropic" component given in lAbdo et al 



(j2009dl ) which investigates the me dium-latitude diffuse emission. We note that the isotropic diffuse 



components obtained here and in lAbdo et al.l (|2009dl ) include the residual background and thus 



should be regarded as an upper limit of the true extragalactic diffuse 7-ray emission. We also note 
that the adopted IC model affects the spectral shape and the intensity of our isotropic component, 
whereas it does not affect the emissivity significantly; modifying the IC emission by ±50 % changes 
the offsets by 6 %-7 %, but alters the slopes less than 3 % except the lowest two energy bands. A 
detailed study of the extragalactic diffuse emission and the residual background using data for a 
larger sky area is underway and will be published elsewhere (A. A. Abdo et al. 2009, in preparation). 
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So far, we have been neglecting the contribution from CR interactions with ionized hydrogen 
(Hll). The low-density ionized gas is unobservable, but can be inferred froni dispe rsion measures 
of pulsar signals in the radio band. According to the model of lCordes fc: Lazid (|2002l ). in the region 
we are studying, A^(HII) is only (1-2) x lO^'^ cm~^ and fairly smooth. We thus conclude that the 
contribution from ionized gas does not affect the obtained emissivity significantly. 



Discussion 



With the approach described in Section 3, we succeeded in decoupling diffuse 7-rays related 
to the local atomic gas from point sources, the IC emission, and the isotropic diffuse compo- 
nent. The derived differential 7-ray emissivity from the local atomic gas is given in Figure 5. 
The systematic uncertainty of the effective area of the response we used (P6_V3_DIFFUSE) 
is estimated to be 10 %, 5 % and 20 % at 100 MeV, 560 MeV, and 10 GeV, respectively, and 
depend on the energy linearly in a logarithmic scale. This systematic uncertainty is compara- 
ble with the statistical error, and is indicated by the shaded area in the figure. The integral 

10^^^ photons s^^ sr^^ H-atom^^ and 
photons ST~^ H-atom~^, respectively, with an additional systematic un- 
^. These values caii be c ompared with those reported by early measurements. 

photons s~^ sr~^ H-atom~^ 



214 emissivity above 100 MeV and 300 MeV is (1.63 ± 0.05) 

215 (0.66 ± 0.02) X 10-26 c-i c.-i tr — -i 

216 certain ty of ^ 10 

217 SAS-2 (IFichtel. Simpson, k Thompsonll 19781 ) obtained about 3 x 10 



218 and COS-B ( Lebrun et al 



19821 ) reported (1.67 ± 0.24) x lO'^^ photons s'^ sr'^ H-at om"^ above 



100 MeV. The EGRET analysis of various directions toward large molecular clouds (IDigel et al. 
2OO1I ) gives (1.65-2.4) xl0~26 photons s"^ sr^^ H-atom^i and (0.71-1.0) x lO'^e photons sr'^ H-atom" 



221 above 100 MeV and 300 MeV, respectively. While most of these early measurements are consistent 

222 with the LAT data, the emissivity obtained by the LAT is much improved in photon statistics and 

223 energy range. 

224 We can give constraints on the local CR spectrum by comparing the obtained emissivity with 

225 the model calculation of interactions between CRs and interstellar matter. Many evaluations of 
the 7-ray production due to CR intera cti ons in the interstellar r nediu r n have been made, in c luding 



Steckerl (Il973l. Il989l). lOermerl (|l986al lbl'l. Isertsch et all (Il993l ). iMoril (|l997l l. IXamae et all ^20od ) 



and iHuang et al.l (|2007l ). In the calculation of neutral pion production and decay 7-rays, most 
authors have computed the 7-ray flux produced through interactions of high-energy CR protons 
with proton targets. The effects of heavy nuclei in both CRs and the target matter are usually 
taken into account as a so-called nuclear enhancement factor (eM) to multiply the proton-proton 
7-ray yield. Although the predicted 7-ray spectra from proton-proton i nteractions assumin g the 



233 same CR proton spectrum agree well (< 10 %) among these works (e.g.. iKamae et al.ll2006l ). the 



234 nuclear enh ancem e nt fac tors differ by u p to ^ 30 %; th e factors range from 1.45 to 1.80- 2.0 as 



235 compi led by Mori ( 2009 ). Among them, Dermer ( 1986a| jbl) gives the lowest eM of 1.45, and Mori 

236 ([20091) gives the highest eM of 1.84 at the CR proton kinetic energy of 10 GeV. His hig her value 

237 of eM is attributed to the adoption of recent CR spectral formulae bv lHonda et al.l (|200J) and the 
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23B inclusion of nuclei heavie r than He in both the interstellar medium and the CR spectra. We thus 
239 regard em by iMoril (|2009l ) as the most reliable. 



240 
241 

242 
243 
244 
245 
246 
247 
24B 

249 
250 
251 
252 
253 
254 
255 
256 
257 
25B 
259 



In calcul a ting t he neutral pion production, we used the proton-proton interaction formalism by 
Kamae et al.l ( 20061 ). They gave parameterized formulae of the vr'^ inclusive cross section and decay 
7-ray spectra for arbitrary proton kinetic energies from 0.488 GeV to 512 TeV. We adopted the 
proton local interstellar spectrum (LIS) from the GALPROP model with 54_5gXva rh7S and caL 



culate d the 7-ray spectrum from nucleon-nucleon interactions using formulae given bv lKamae et al 



(j200g) und er the assumption of the nuclear enhancement factor to be 1.84 as a representative value 
of those by iMoril (|2009l ) . In order to calculate the electron brem sstrahlung, we fully u tilized GAL- 
P ROP which calculates the 7-ray spectru m using a formalism bv lKoch &: Motzj (|l959l ) as explained 
by IStrong. Moskalenko. Reimed (|200d ) . 



The predicted emissivity at the solar sytem (Galactocentric radius R = 8.5 kpc and the height 
from the Galactic plane z = kpc) is compared with our LAT measurement in Figure 5, and 
the LIS of proton, electron and positron used in this model calculation is presented in Figure 6. 
Also shown is a compilation of some measurements of proton and electron spectra at the Earth. 
The proton and electron model spectra follow the observed ones above a few tens of GeV; below 
this energy the solar wind lowers the observed fluxes. We note that th e CR electron spe ctrum 
measured by the LAT is somewhat harder than the GALPROP model (lAbdo et al.ll2009d ). but 
the effect on our analysis is negligible. We also note that the bremsstrahlung at around 100 MeV 
has compar able contributions f rom both primary electrons and secondary electrons/positrons, as 
discussed by lPorter et al.l ([20081). Their contributions are included in the computed spectrum shown 
in Figure 5. 



260 The true LIS is somewhat uncertain d ue to solar activity . To model this solar modulation 

261 effect on the CR spectrum, the formula by iGleeson &: AxfordI (jl968l ) generally has been used, in 

262 which a single parameter <j) is introduced. The proton LIS we adopted reproduces the observations 

263 at the Earth with <j) = 450 MV as shown in Figure 6. The same data can also be r eproduced by 

264 a diffe rent formula of the LIS (dotted blue line) with cp = 600 MV as described by IShikaze et al 



265 (120071 ). We thus regard the difference between two models as representing the uncertainty of the 

266 LIS; they agree well (< 10 %) above 10 GeV and differ by ~ 20 % at about 3 GeV. This affects the 

267 calculated emissivity by 20 % and < 10 % at about 100 MeV and above 1 GeV, respectively (e.g.. 



Morilll997l ). We note that the LIS we adopted better reproduces the observed proton spectrum 



269 above 20 GeV. Although the true LIS below 1 GeV is highly uncertain, these CRs do not contribute 

270 to the 7-ray emissivity above 100 MeV significantly. 

271 As shown in Figure 5, the emissivity measured by the LAT agrees with the prediction from 

272 the assumed LIS and the recent estimate of en at the 10 % level, which is comparable with the 

273 statistical error and the current systematic uncertainty of the LAT response. For re ference, we 

274 also show the 7-ray emissivity model obtained with the lowest eM (among references in lMorill2009l ) 

275 of 1.45 that gives a predicted emissivity lower than that observed in 100 MeV-10 GeV. Since the 
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276 nucleon-nucleon component is dominant in the emissivity spectrum especially above 1 GeV, the 

277 observed agreement between the LAT data and the model calculation (with the latest estimate of 

278 eM) indicates that CR nuclei in the vicinity of the solar system in regions observed have spectral 

279 distributions and intensities close to those of the LIS inferred from measurements at the Earth 

280 within ~ 10 %. Although the constraint is rather weak, the agreement down to 100 MeV also 

281 suggests that highly uncertain low-energy (below a few hundred MeV) CR electron and positron 

282 spectra are compatible with our assumption (GALPROP model with 54_5gXvarh7S) shown in 

283 Figure 6. 

284 5. Summary and Conclusions 

285 Wc report the observation of diffuse 7-rays in a mid-latitude region in the third quadrant using 

286 data from the first six months of Fermi LAT science observations. The region is away from the 

287 Galactic plane and the Galactic center, and contains no known large molecular cloud. Most of the 

288 atomic hydrogen is within 1 kpc of the solar system, and thus the region is suitable for studying the 

289 7-ray emissivity of the local atomic gas and CR spectra in the neighborhood of the solar system. 

290 Thanks to the excellent performance of the LAT and recent developments of the CR propagation 

291 code and the interstellar radiation field model in GALPROP, we reliably estimated and subtracted 

292 the contribution from point sources and IC emission. The residual 7-ray intensities exhibit a 

293 linear relationship with the atomic gas column densities from 100 MeV to 10 GeV, indicating that 

294 non-isotropic 7-rays are produced through interactions of CRs with interstellar atomic gas. The 

295 measurement of the emissivity of local atomic hydrogen has already surpassed those by past missions 

296 in photon statistics and the energy range. It agrees with the prediction from CR spectra assumed, 

297 indicating that the CR nuclei spectra in the vicinity of the solar system in regions analyzed are 

298 close to the LIS inferred from direct measurements at the Earth within ~ 10 %. Low energy CR 

299 electron/positron spectra are suggested to be compatible with our assumption. 
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312 Some of the results in this paper have been derived using the HEALPix (|G6rski et al.l |2005| ) 

313 package. 
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Fig. 1. — Fermi LAT 7-ray count maps {E > 100 MeV) of regions we analyzed. Maps are in 
Cartesian projection with 0.5° x 0.5° binning. The north and the south regions are shown in the 
upper panel and the lower panel, respectively. In these regions there are 52 sources in the LAT 6 
month sourc e list, of which 13 s ources (indicated by green crosses) are included in the LAT Bright 
Source List (|Abdo et ahlboOObl ^. 
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Fig. 2. — Diffuse 7-ray spectrum in regions analyzed {I from 200° to 260° and \h\ from 22° to 60° ) 
after masking point sources witli circular regions of 1° radius. IC emission predicted by GALPROP 
(54_5gXvarh7S) and the residual point source contributions estimated from spectral parameters 
given in the LAT 6 month source list are shown by dotted and solid histograms, respectively. The 
horizontal and vertical error bars indicate the energy ranges and 1 a statistical errors, respectively. 
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Fig. 3. — Column density maps of atomic hydrogen derived from the LAB survey under the 
assumption of a uniform spin temperature of 125 K. 
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Fig. 4. — Correlation of the (IC and point-sources subtracted) 7-ray intensities and the HI column 
densities in four representative bands. The map of A^(HI) (Figure 3) is convolved with the LAT 
PSF of the corresponding energy range. The horizontal and vertical error bars indicate the ranges 
of the column density and the 1 a statistical errors, respectively. Data in high energy range (above 
1 GeV) are rebinned to have more than 10 7-ray counts in each bin. 
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Fig. 5. — Differential 7-ray emissivity from the local atomic hydrogen gas compared with the 
calculated 7-ray production. The horizontal and vertical error bars indicate the energy ranges and 
1 a statistical errors, respectively. The assumed interstellar proton, electron and positron spectra 
are shown in Figure 6. Estimated systematic errors of the LAT data are indicated by the shaded 
area. A nucleus enhancement factor eM of 1.84 is assumed for the calculation of the 7-rays from 
nucleon-nucleon interactions. Dotted lines indicate the emissivities for the case of cm = 1-45, the 
lowest values in the referenced literature. 
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Fig. 6. — LIS of protons (thick blue line), electrons (thick red line) and positrons (thin red line) 
used to calculate the 7-rav emis s ivity ra o del in Figure 5 with a cor npilation of direct me asurements 
at the Earth; see 



Alcaraz et all (12000^) 



Sanuki et al 



_ (jgOOm and IShikaze et alJ (|2007l l for proton 

data, and lAlcaraz et alJ (j2000bl ) and Bar wick et al. ( 19981 ) for electron dat a. Modulated s pectr um 
with (j) = 450 MV is given by thin blue line. The proton LIS adopted by IShikaze et al.l (|2007l ) is 
shown by dotted blue line. 
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Table 1. A summary of fit parameters with 1 Sigma statistical errors. 



Energy 
GeV 


Offset 

10"'' MeV^ s"^ cni-2 sr-^ MeV"^ 


Slope 

10-2* MeV^ s"^ sr"^ MeV"^ 


xVdof^ 


Data/Model counts 


0.10-0.14 


15.40 ±0.54 


1.04 ±0.14 


21.90/13 


11799/11678 


0.14-0.20 


17.10 ± 0.40 


1.67 ±0.10 


18.12/14 


27891/27738 


0.20-0.28 


16.70 ± 0.36 


1.91 ±0.09 


13.47/14 


31718/31564 


0.28-0.40 


15.83 ±0.36 


2.11 ± 0.10 


16.92/14 


28987/28850 


0.40-0.56 


13.81 ± 0.39 


2.33 ± 0.10 


6.65/14 


22718/22073 


0.56-0.80 


12.57 ±0.41 


2.20 ±0.11 


16.17/15 


16137/16063 


0.80-1.13 


11.44 ±0.44 


2.17 ± 0.12 


12.26/15 


11421/11368 


1.13-1.60 


10.23 ± 0.49 


1.88 ±0.13 


9.06/14 


7364/7327 


1.60-2.26 


9.25 ± 0.54 


1.72 ±0.15 


13.16/14 


4765/4733 


2.26-3.20 


8.44 ± 0.58 


1.15 ± 0.16 


10.56/12 


2764/2747 


3.20-4.53 


7.12 ± 0.64 


1.10 ±0.17 


15.07/11 


1733/1712 


4.53-6.40 


6.44 ± 0.75 


1.12 ±0.21 


11.35/11 


1158/1144 


6.40-9.05 


5.51 ±0.77 


0.71 ±0.21 


9.92/9 


678/664 



'^Degree of freedom 

''Data and model total counts after masking point sources with circular regions of 1° radius. We believe that the 
small (< 1 %) excesses of the data counts over the model counts are due to unresolved point sources or interstellar 
matter not traced by 21 cm line surveys. 



